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Abstract: Till the middle of 2004 it appeared that the unitarity relation + + 

= 1 might not hold at the 2.3a level. At that time, however, \Vus\ was inferred 
from old experimental data. Since then, a large experimental and theoretical effort has 
been invested leading to a removal of the problem. Thanks to the new and improved 
measurements by BNL-E865, KLOE, KTeV, ISTRA+ and NA48, the K^s decay rate moved 
up so that \Vus\ is now consistent with unitarity. On the theory side, much progress has 
been made in order to tame the systematic uncertainties related to the computation of the 
form factors. 

This joint progress allowed to assess the validity of the CKM unitarity relation at the level 
of less than 1%. Recent measurements of kaon decays contributing to the determination 
of \ Vus\ \Vus\/\Vud\ are summarized, and up-to-date evaluations of |14s|/-|-(0) and \ Vus\ are 
presented. In addition, we discuss the sensitivity of Kg^ and K£2 decays to various scenarios 
of physics beyond Standard Model. 
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1. Introduction 

I report on precise tests of the Standard Model (SM) with kaon decays using world data. 
From the experimental information on down- to up-quark transitions (such as d ^ u, 
s ^ u and b ^ u), we access the effective dimesion-six operators of the form, DTiUiT2V, 
with D (U) being a generic ^^down" {^''up") flavor, and i = e, ^, r. Their effective cou- 
pling are parametrized as the SM contribution G'j^ \ Vud\'^, plus a possible new physics 
terms, G'^pe^p. Since the dimension-six operators are not protected by gauge invariance 
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the possible effects of non-decouphng are proportional to (1 + M^/A^p). The effects 
of these non-standard contributions cannot be very large, but are possibly detectable in 
high-precision experiments. 

A convenient strategy to measure these effects against the SM parameters, G'p and 
iVf/ul, rely on the Cabibbo universality hypothesis (or unitarity constraint): 

GIkm = Gl (or iKdl' + \Vus\^ + iKbl' = 1 and Gp = G^) , (1.1) 

where Gl^^ = + |K«P + iKfeP) , and G^ = 1.166371(6) x 10-5GeV-^ as 

extracted from the accurate measurement of the muon lifetime [||. 

I report on the progress related to the verification of the unitarity relation ( |1 . Ij) . As 
we shall see the current accuracy of the CKM unitarity relation ( |1.1D , is at the 0.1% level, 
becoming an important constraint to the model builders of scenarios beyond SM physics. 
Kii2, and K^2 decays offer possibly the cleanest way to test us transitions, ud transitions 
are precisely measured in superallowed nuclear /3-decays. The most recent determination 
of Kd is |Kd|=0.97418± 0.00026, g. 

This report is organized as follows. The phenomenological farmework needed to de- 



scribe Ki'i and K^2 decays is briefly recapitulated in Section |1 . 1| . Section y is dedicated to 
the combination of the experimental data. The results and interpretations are presented 
in Section |3|. 

1.1 Kt2 and K^^ phenomenology 

For {T^n) amplitudes, we introduce the following QCD parameters 
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mi 



(0|s7^75u|ir+ {p)) = iV2fKP^, (0|s75n|K+) = -iV2fK , (1.2) 

rus + ruu 

For amplitudes, we define the following form factors 



V2 



(1.3) 



where t = {p — k) and 



/Vf2 _ jur2 

(vr+ {k) \-su\K' (p)) = ^—h (t) (1.4) 

V2 [ms - niu) 

The SM gives the following relations for the K^^ and K12 decay rates: 

r(i^^3(7)) = ifefC-^^- \Vus\'f+{Of 4(A+,o) (1 + 5ff,(2) + 5em') , (1-5) 



± ] 



r(^^2(7)) 



Vus 



V, 



ud 



K 



x(l + 5cm) , (1.6) 



where Ck = 1 (1/2) for the neutral (charged) kaon decay. /^(A+^o) is the phase space 
integral which also includes the form factors parameterized by A+ o- The universal short- 
distance electromagnetic correction. Sew = 1.0232(3), has been computed at n = Mp in 



2 



ref. [P, while the long-distance electromagnetic corrections, Sem = 0.9930(35) P and (5^, as 
well as the isospin-breaking ones, (^^jy^j)' have been recently revised in ref. ||6| (see table |l]). 









K% 




2.36(22) 


2.36(22) 


+0.57(15 
+0.08(15) 
+0.80(15) 
+0.05(15) 



Table 1: Summary of the isospin-breaking factors [Q, ^ 
with correlations for 5^(%): 



/ 1. 0.11 0.78 -0.12 \ 
1. -0.12 0.78 
1. 0.11 
1. 



(1.7) 



The remaining quantities, /-i-(0), the vector form factor at zero momentum transfer [q^ = 
{pk — Pn)'^ = 0], and fK/f-K, the ratio of the kaon and pion decay constants contains 
the non-perturbative QCD information on the flavor SU(3) breaking effects arising in the 
relevant hadronic matrix element. For more details see Sec. |3.2.1 . 



1.2 Form Factors Parameterizations: Lattice QCD, ChPt and Dispersion rela- 
tions 

To determine Vus , we have to determine the integral over phase space of the Dalitz density, 
which depends on the form of f+^o{t), the form factors (FF) in eq. ( |1.3D . To reduce 
uncertainties, it would be convenient to have theoretical information on their t dependence. 
Our theorical knowledge is still poor. ChPt and Lattice QCD could be very useful but 
predictivity of ChPt is limited by our knowledge of Low Energy Constants(LECs), while 
lattice calculations still have large uncertainty. 

For this reason, present estimates of the phase space integral mainly rely on measure- 
ments. Li the physical region, it is a good approximation to use a quadratic t dependence 
of the FFs such as 



f+At) 



Note that t = {px — PwY 



f+At) 

/+,o(0) 



l + A' 



,0 ^2 



+,0 



(1.8) 



mj^ + — 2m;^£'^, therefore the FFs depends only on 
Ej^. The FF parameters can thus be obtained from a fit to the pion spectrum which is of 
the form g^E^^) x /(E^^)'^. Unfortunately t is maximum for £^^=0, where g{ET^) vanishes. 
Still, experimental information about /+ are quite accurate and a pole parametrization, 
A" ~ 2(A'^)^, looks confirmed from present data (see later). Theory would also support this 
hypothesis. The situation of /o instead is completely open. The main problem is that Aq 
cannot be detected from the data and we can not discriminate among different assumptions 
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such as hnear, pole or quadratic fits. In turn, these model ambiguities induce a systematics 
uncertainty for |V^s| even though data for partial rates by itself are very accurate. For 
this reason, hints from theory are welcome and it is advisable to use model-independent 
parametrisations as possible. 

The vector and scalar form factors f-^-/Q{t) in eq. ( |1.3| ) are analytic functions in the 
complex t~plane, except for a cut along the positive real axis, starting at the first physical 
threshold tthi = {Mk + Af^r)^, where their imaginary parts develop discontinuities. They 
are real for s < tthi- 

Cauchy's theorem implies that f^/oit) can be written as a dispersive integral along the 
physical cut 

oo ^ 

f(t) = - f ds' , j'^f^^ ^ + subtractions . (1.9) 
TT J [s' —t — iO) 

hh 1 

where all possible on-shell intermediate states contribute. A number of subtractions are 
needed to make the integral convergent. Particularly appealing is an improved dispersion 
relation recently proposed in ref. [^, which reads 

/.(*) = ( UtcT) e(*-')'"") * , mt) = & r ^ , , t'r' , ,, (1.10) 

V / Jt^^ X {x- tcT) [x-t- le) 

Here 0(x), the phase of foit), is taken from Kir scattering. This dispersive form has been 
solved in terms of the subtractions at t = (where by default, /o(0) = 1 is known) and at 
tcT = ("T-zf — "^^)) where the QCD soft-pion theorem, known as Callan-Treiman relation, 
implies 

Mt ^ (rn\ - ml)) = + (3-5 ± 8.0 lO'^g) (1.11) 

By exploiting this form on the spectrum of Ki^ data, we have to estimate one unknown 
parameter {fo{t = {rriK — rriT^)'^) or fK/fw/f+iO)), whereas, eq. (|1.10 ) within its theoretical 



uncertainty non-trivially contraints the coefficients of the Taylor's expansion in eq. ( |1.8| )(for 
axample for the ffi'st two derivatives, we have from ref. P Ag = (Ag)^ — 2m^/tcTG'{0) = 
(Ao)^ + (4.16±0.50) X 10~^). By the experiemtal information of fo{t = (mK — m.,,)'^), Callan- 
Treiman theorem allows for a consistency checks with lattice QCD, where fK/fn/f+{0) can 



be estimated (see section 3. 2. 3) 



2. Data Analysis 

We perform fits to world data on the BRs and lifetimes for the Kl and K^, with the 
constraint that BRs add to unity [0. This is the most satisfactory way of making use of 
the new measurements. 

2.1 Kl leading branching ratios and tl 

Numerous measurements of the principal Kl BRs, or of various ratios of these BRs, have 
been published recently. For the purposes of evaluating |T4s|/+(0), these data can be used 
in a PDG-like fit to the Kl BRs and lifetime, so all such measurements are of interest. 
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KTeV has measured five ratios of tlie six main Kl BRs 0. The six channels 
involved account for more than 99.9% of the Kl width and KTeV combines the five 
measured ratios to extract the six BRs. We use the five measured ratios in our anal- 
ysis: BR(K^3/Ke3) = 0.6640(26), BR(7r+7r^7r°/i^e3) = 0.3078(18), BR(7r+7r^//Ce3) = 
0.004856(28), BR{^n^/Ke3) = 0.4782(55), and BR(27rO/37rO) = 0.004446(25). The errors 
on these measurements are correlated; this is taken into account in our fit. 

NA48 has measured the ratio of the BR for Kez decays to the sum of BRs for all 
decays to two tracks, giving BR(Ke3)/(l - BR(37r°)) = 0.4978(35) From a separate 

measurement of BR(Kl ^ 37r°)/BR(K5 27r°), NA48 obtains BR(37rO)/rL = 3.795(58) 

/xs-i m. 

Using (j) KlKs decays in which the Ks decays to 7r"'"7r~, providing normalization, 
KLOE has directly measured the BRs for the four main Kl decay channels |14]. The 
errors on the KLOE BR values are dominated by the uncertainty on the Kl lifetime tl', 
since the dependence of the geometrical efficiency on tl is known, KLOE can solve for tl 
by imposing BR(i^'i — > x) = 1 (using previous averages for the minor BRs), thereby 
greatly reducing the uncertainties on the BR values obtained. Our fit makes use of the 
KLOE BR values before apphcation of this constraint: BR(i^e3) = 0.4049(21), BR(E:^3) = 
0.2726(16), BR(/Ce3) = 0.2018(24), and BR(Ke3) = 0.1276(15). The dependence of these 
values on tl and the correlations between the errors are taken into account. KLOE has 
also measured tl directly, by fitting the proper decay time distribution for Kl 37r^ 
events, for which the reconstruction efficiency is high and uniform over a fiducial volume 
of ~0.4Al. They obtain tl = 50.92(30) ns |l|. 

There are also two recent measurements of BR(7r"'"7r~/i^£3), in addition to the KTeV 
measurement of BR(7r^7r~/i^^e3) discussed above. KLOE obtains BR(7r~''7r~/ir^3) = 
7.275(68) X 10-3 [||], while NA48 obtains BR(7r+7r-/iCe3) = 4.826(27) x 10"^ @. All 
measurements are fully inclusive of inner bremsstrahlung. The KLOE measurement is fully 
inclusive of the direct-emission (DE) component, DE contributes negligibly to the KTeV 
measurement, and a residual DE contribution of 0.19% has been subtracted from the NA48 
value to obtain the number quoted above. For consistency, in our fit, a DE contribution 
of 1.52(7)% is added to the KTeV and NA48 values. Our fit result for BR(7r+7r-) is then 
understood to be DE inclusive. 

In addition to the 14 recent measurements listed above, our fit for the seven largest 
Kl BRs and lifetime uses four of the remaining five inputs to the 2006 PDG fit and the 
constraint that the seven BRs add to unity. The results are given in Table ^. 

The evolution of the average values of the BRs for Kliz decays and for the important 
normalization channels is shown in Fig. ||. 

Our fit gives xV^df = 20.2/11 (4.3%), while the 2006 PDG fit gives xV^df = 14.8/10 
(14.0%). The differences between the output values from our fit and the 2006 PDG fit are 
minor. The poorer value of x^/ndf for our fit can be traced to contrast between the KLOE 
value for BR(37r'') and the other inputs involving BR(37r'') and BR(7r*'7r'') — in particular, 
the PDG ETAFIT value for BR(7r07r°/7r+7r-). The treatment of the correlated KTeV and 
KLOE measurements in the 2006 PDG fit gives rise to large scale factors for BR(ife3) and 
BR(37r'^); in our fit, the scale factors are more uniform. As a result, our value for BR(i^e3) 
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I—' O V O TTI mT /^V 


Value 


I? 
o 




U.4lU0Do^ ( 4j 


i.i 




n 97n47/'71 


1 1 
1 . 1 


BRf37r°l 


iq5r)7('86') 


1.2 


BR(7r+7r-7rO) 


0.12542(57) 


1.1 


BR(7r+7r-) 


1.9966(67) X 10-3 


1.1 


BR(27rO) 


8.644(42) X 10-^ 


1.3 


BR(77) 


5.470(40) X 10"^ 


1.1 


TL 


51.173(200) ns 


1.1 



Table 2: Results of fit to Ki BRs and lifetime 
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Figure 1: Evolution of average values for main Kl BRs. 

has a significantly smaller uncertainty than does the 2006 PDG value. 
2.2 Ks leading branching ratios and ts 



KLOE has published |17| a measurement of BR(Er5 — > neu) that is precise enough to 



contribute meaningfully to the evaluation of |T4s|/+(0). The quantity directly measured 
is BR(7rez/)/BR(7r+7r-). Together with the published KLOE value BR(7r+7r-)/BR(7r°7rO) 
= 2.2459(54), the constraint that the Ks BRs must add to unity, and the assumption of 
universal lepton couplings, this completely determines the Ks BRs for tt+tt-, vr'^vr'', K^s, 
and Kf,3 decays Q. In particular, BR{Ks ireu) = 7.046(91) x 10"^. 

NA48 has recently measured the ratio T{Ks ■Kev)/T{KL -nev) = 0.993(26) (22) 



1 19]. The best way to include this measurement in our analysis would be via a combined 
fit to Ks and K^ branching ratio and lifetime measurements. Indeed, such a fit would be 
useful in properly accounting for correlations between Ks and K^ modes introduced with 
the preliminary NA48 measurement of T{Kl Stt*^), and more importantly, via the PDG 
ETAFIT result, which we use in the fit to Kl branching ratios. At the moment, however, 
we fit Ks and Kl data separately. NA48 quotes BR(i^5 -jrev) = 7.046(180) (160) x lO"'^; 
averaging this with the KLOE result gives BR{Ks — > vrei^) = 7.046(84) x 10~^, improving 
the accuracy on this BR by about 10%. 

For TKg we use 0.8958 x 10"^*^ s, where this is the non-CPT constrained fit value from 
the PDG, and is dominated by the 2002 NA48 and 2003 KTeV measurements. 
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I—' O T* O TTl JTT r^t^ 




o 




Do.ODy\^iio ) /O 


1.1 




20 644('80')% 


1.1 


BR(7r7r7r) 


5.5953(308)% 


1.0 


BR(i^e3) 


5.0780(258)% 


1.2 


BR(i^;.3) 


3.3650(271)% 


1.7 


BR(7r7r07r°) 


1.7495(261)% 


1.1 




12.3840(193) ns 


1.7 



Table 3: Results of fit to BRs and lifetime 



2.3 leading branching ratios and 

There are several new results providing information on rates. These results are mostly 
preliminary and have not been included in previous averages. 

NA48/2 has recently published measurements of the three ratios BR(Ere3/7r7r'^), BR(i^^3/7r7r'^ 
and BR{K^3/ Kes) [20|. These measurements are not independent; in our fit, we use the 
values BR(i^e3/vr7r°) = 0.2470(10) and BR{K^s/7nT^) = 0.1637(7) and take their correla- 
tion into account. ISTRA+ has also updated its preliminary value for BK{Ke3/7nT^). They 
now quote BR{Ke3/inT^) = 0.2449(16)113]. 

KLOE has measured the absolute BRs for the i^es and -ftT^s decays and a very 
precise measuremet of BR(i^^2)||Hl- In (/> — > K^K~ events, decays into fiu or ttt:^ 
are used to tag a K~ beam, and vice versa. KLOE performs four separate measurements 
for each i^£3 BR, corresponding to the different combinations of kaon charge and tagging 
decay. The final averages are BR(ire3) = 4.965(53)% and BR(ir^3) = 3.233(39)%. 

Very recently KLOE has also measured the absolute branching ratio for the vrvr'' decay 
with 0.5% accuracy. The KLOE preliminary result, is BR(7r7r°) = 0.20658(112) |2|]. 

Our fit takes into account the correlation between these values, as well as their depen- 
dence on the lifetime. The world average value for t± is nominally quite precise; the 
2006 PDG quotes r± = 12.385(25) ns. However, the error is scaled by 2.1; the confidence 
level for the average is 0.17%. It is important to confirm the value of t±. The two new 
measurements from KLOE, t± = 12.364(31)(31) ns and t± = 12.337(30)(20) ns[||] with 
correlation 34%, agree with the PDG average. 

Our fit for the six largest BRs and lifetime makes use of the results cited above, 
plus the data used in the 2006 PDG fit, except for the Chiang '72 measurements for a total 
of 26 measurements. The six BRs are constrained to add to unity. The results are shown 
in Table |. 

The fit quality is poor, with x^/ndf = 42/20 (0.31%). However, when the five older 
measurements of t± are replaced by their PDG average with scaled error, x'^/ndf improves 
to 24.3/16 (8.4%), with no significant changes in the results. 

Both the significant evolution of the average values of the -ftr^3 BRs and the effect of 
the correlations with BR(7r7r'^) are evident in Fig. |2[ 
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0.62 



0.63 
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0.035 
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-■- 
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0.64 0.195 



0.205 



0.215 



Figure 2: Evolution of average values for main BRs. 



2.4 Measurement of BR(i^e2)/BR(K^.2) 



Experimental knowledge of Kf>2/ K^2 has been poor so far. The current world average of 
Rk = BR(Ke2)/BR(i^^2) = (2.45 ± 0.11) x 10"^ dates back to three experiments of the 
1970s [^] and has a precision of about 5%. The three new preliminary measurements were 
reported by NA48/2 and KLOE (see Tab. |): A prehminary result of NA48/2, based on 
about 4000 events from the 2003 data set, was presented in 2005 [^7|. Another pre- 
liminary result, based on also about 4000 events, recorded in a minimum bias run period 
in 2004, was shown at KAONOT p^ . Both results have independent statistics and are also 
independent in the systematic uncertainties, as the systematics are either of statistical na- 
ture (as e.g. trigger efficiencies) or determined in an independent way. Another preliminary 
result, based on about 8000 K(.2 events, was presented at KAON07 by the KLOE collab- 
oration Both, the KLOE and the NA48/2 measurements are inclusive with respect 
to the final state radiation bremsstrahlung contribution. The small contribution of Ki2^ 
events from direct photon emission from the decay vertex was subtracted by each of the 
experiments. Combining these new results with the current PDG value yields a current 
world average of 

Rk = (2.457 ± 0.032) x 10"^ (2.1) 



in very good agreement with the SM expectation and, with a relative error of 1.3%, a factor 
three more precise than the previous world average. 
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Rk [10-5] 


PDG 2006 H] 


2.45 ±0.11 


NA48/2 prel. ('03) |7 


1 2.416 ± 0.043 ± 0.024 


NA48/2 prel. ('04) ^ 


I 2.455 ± 0.045 ± 0.041 


KLOE prel. [|§ 


2.55 ±0.05 ±0.05 


SM prediction 


2.477 ±0.001 



Table 4: Results and prediction for Rk = BR{Ke2)/^R{K^2)- 



2.5 Measurements of K^s slopes 

2.5.1 Vector form factor slopes from K^^ 

For i^e3 decays, recent measurements of the quadratic slope parameters of the vector form 



factor (AV, A'l) are available from KTeV |g], KLOE ISTRA± ||, and NA48 ||. 

We show the results of a fit to the Kl and data in the first column of Table and 
only the Kl data in the second column. With correlations correctly taken into account, 
both fits give good values of x^/ndf. The significance of the quadratic term is 4.2(T from 
the fit to all data, and 3.5a" from the fit to Kj^ data only. 





Kl and K data 


Kl data only 




4 measurements 


3 measurements 




xVndf = 5.3/6 (51%) 


xVndf = 4.7/4 (32%) 


A'+ X 10^ 


25.15 ±0.87 


24.90 ± 1.13 


A'l X 103 


1.57 ±0.38 


1.62 ±0.46 


P(AV,A±") 


-0.941 


-0.951 


IiK%) 


0.154651(236) 


0.154560(307) 




0.159005(241) 


0.158912(315) 



Table 5: Average of quadratic fit results for Ke^ slopes 



Including or excluding the K slopes has little impact on the values of A^ and X'!^; in 
particular, the values of the phase-space integrals change by just 0.07%. 

KLOE, KTeV, and NA48 also quote the values shown in Table |6| for My from pole 
fits to Kl e3 data. The average value of My from all three experiments is My = 875.3 ± 
5.4 MeV with x^/ndf = 1.80/2. The three values are quite compatible with each other and 
reasonably close to the known value of the K^*{SQ2) mass (891.66 ±0.26 MeV). The values 
for A^ and A" from expansion of the pole parameterization are qualitatively in agreement 
with the average of the quadratic fit results. More importantly, for the evaluation of the 
phase-space integrals, using the average of quadratic or pole fit results gives values of /(iCgg) 
that differ by just 0.03%. No additional error needs be assigned to account for differences 
obtained with quadratic and pole parameterizations for the vector form-factor slope. 
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Experiment 


Mv (MeV) 


(My) = 875.3 ± 5.4 MeV 


KLOE 


870 ± 6 ± 7 


xVndf = 1.80/2 


KTeV 


881.03 ±7.11 


A'^ X 10^ = 25.42(31) 


NA48 


859 ± 18 


A'l = 2 X 
I{K^^) = 0.154695(192) 



Table 6: Pole fit results for K^^ slopes 



2.5.2 Scalar and Vector form factor slopes from K^s 



For JC^3 decays, recent measurements of the slope parameters (A'_,_, A'_^, Aq) are available 
from KTeV ||], KLOE ||], 1STRA+ |3|, and NA48 [|6|. Note that it is not possible, 



because of correlations, deduce the presence of quadratic term in /o(i) from the decay 
spectra. For the same reason fits with a linear parametrization give a wrong result for the 
slope Aq. 

Figure ^ shows the 1-a contours from all the experimantal results (i^es and K^^). It is 
immediately clear from the figure that the new NA48 results are difficult to accommodate. 
Performing the combinaton with and without the NA48 results for the iT^a form-factor 
slopes included we obtain fit probability values of 1 x 10~^ and 22.3% respectively. The 
results of the combination are listed in Table |^. 
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KTeV 



ISTRA+ 



NA48 




25 - 



I'+xlO' 




" I'+xlO^ 




: c 

1 1 1 1 1 1 


1 

1 , . , 



10 



15 



AqXIO" 



Figure 3: l-cr contours for A^, A" , Aq determination from lSTRA+(pink ellipse), KLOE (blue 
ellipse), KTeV(red ellipse), NA48(green ellipse), and world average with(filled yellow ellipse) and 
without (filled cyan ellipse) the NA48 i^T^s result. 



The value of x^/ndf for all measurements is terrible; we are forced to quote the results 
with scaled errors. This leads to errors on the phase-space integrals that are ~60% larger 
after inclusion of the new Kf^^ NA48 data. 
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Kl and K 




i^L only 




Measurements 


16 




11 




X /ndf 


54/13 (7 X 10 ' 


) 


33/8 (8 X 10 ^) 


1 


X lO'^ 


24.920 ± 1.105 (5 = 


: 1.4) 


24.011 ± 1.544 (S = 


: 1.5) 


A" X 10^ 


1.612 lb 0.447 (6 = 


1.3) 


1.974 lb 0.622 (5 = 


1.6) 


Ao X 10^ 


13.438 ± 1.19 (6 = 


1.9) 


11.682 lb 1.238 (6 = 


: 1.7) 


/ \ f \ n \ 

P(A+,A+) 


—0.944 




—0.966 




/5(A+, Ao) 


+U.3zo 




+0.715 




p(A+, Aoj 


n /I Qo 
— U.4oy 




— U.Dyo 






0.154566(292) 




0.154354(389) 






0.158918(300) 




0.158700(400) 






0.102123(312) 




0.101643(424) 






0.105073(321) 




0.104578(437) 






+0.63 




+0.89 





Table 7: Averages of quadratic fit results for i^es and slopes. 



We have checked to see if the NA48 data might show good consistency with the 
results of some other experiment in a less inclusive average. Fitting to only the results 
from KTeV, NA48, and ISTRA+ gives xV^df = 27.5/6 (0.01%). Fitting to only the K^^z 
results from KTeV, NA48 gives xV^df = 11.6/3 (0.89%). The consistency of the NA48 
data with these other measurements appears to be poor in any case. 

The evaluations of the phase-space integrals for all four modes are listed in each case. 
Correlations are fully accounted for, both in the fits and in the evaluation of the integrals. 
The correlation matrices for the integrals are of the form 

+1 +1 p p 
+1 +1 p p 

P p +1 +1 

p p +1 +1 

where the order of the rows and columns is -fiT^g, iT^, i^^^g, -^^^3, and p = p{Ie3,I^i3) as 
listed in the table. 

Adding the K^^ data to the fit does not cause drastic changes to the values of the 
phase-space integrals for the i^es modes: the values for I{K^^) and I{K^) in Table |^ are 
qualitatively in agreement with those in Table ^. As in the case of the fits to the Ke3 data 
only, the significance of the quadratic term in the vector form factor is strong (3. 60" from 
the fit to all data). 

3. Physics Results 

3.1 Determination of /+(0)Ks and Ks/Kd x /t^//^ 

This section describe the results that are independent on the theoretical parameters /+(0) 
and fK/fn- 
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3.1.1 Determination of /+(0)T4 



The value of /+(0)Kjs has been determined from L5 using the world average values reported 
in |2| for lifetime, branching ratios and phase space integrals, and the radiative and SU(2) 



breaking corrections discussed in section 1.1. 

The results are given in Table |, and are shown in Fig. |^ for Kl 

— > vr^i/, and for the combination. 



mode 


/+(o)v;. 


% err 


BR 


r 


A 


Int 


Kl^ 




0.21625(60) 


0.28 


0.09 


0.19 


0.15 


0.09 


Kl^ 




0.21675(66) 


0.31 


0.10 


0.18 


0.15 


0.15 


Ks^ 


Txev 


0.21542(134) 


0.67 


0.65 


0.03 


0.15 


0.09 


K^- 


> Tieu 


0.21728(84) 


0.39 


0.26 


0.09 


0.26 


0.09 






0.21758(111) 


0.51 


0.40 


0.09 


0.26 


0.15 


averag 


56 


0.21661(46) 













Table 8: Summary of f+{0)Vus determination from all channels. 



The average, |Ks|/+(0) = 0.21661(46), has 
an uncertainty of about of 0.2%. The results from 
the five modes are in good agreement, the fit prob- 
ability is 58%. In particular, comparing the values 
of /+(0)Kis obtained from K^^ and Kf^ we obtain 
a value of the SU(2) breaking correction 



0.215 



0.2175 



lavi 

'net 



Kaon WG 



2.86(39)% 



in agreement with the CHPT calculation reported 
in table |l| (5f^(2) = 2.36(22)%. 



I 



f+(0)xVus 



K 



Le3 



K 



K 



Se3 



e3 



fi3 



0.215 



0.2175 



Figure 4: Display of /+(0)Kis for all 
channels. 



3.1.2 Determination of Vus/Vud x Ik/ fw 

Another determination of \Vus\ is obtained from 

decays. The most important mode is — > 

which has been recently updated by KLOE, 
so that the relative uncertainty is now about 0.3%. 
To reduce hadronic uncertainties, in eq. ( |1.6| ) we have introduced the ratio T{K~^ — > 
//+z^)/r(7r+ ^ fi+iy). 

Using the world average values of BR(i^^ Z^^'^) ai^d of given in section|2| and the 
value of r(7r^ — > fi^v) = 38.408(7)/us-i from |6[ we obtain: 

Vus/Vud X fx/fn = 0.27599 ± 0.00059 

3.2 The parameters /+ and Jk/Itv 

For the time being, such a highly precise measurement could not be translated to a similar 
error on the \ Vus\ determination and therefore to the test of CKM unitarity. The obstacle 
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is obviously the difficulty to keep the theoretical uncertainties in /+(0) and fK/f-w at the 
per-mil level. 



3.2.1 /+(0) determination 

In eq. ( |l.5D , /+(0) is defined in the absence of electromagnetic corrections and I-spin 
breaking effects. While QCD is flavor blind the mass differences ^ rrid ^ mg results 
in /+(0) being different from unity and also being different for charged and neutral kaons. 
In the following, by common convention, /+(0) is refers to neutral kaons, the difference for 
charged kaons is accounted for in the term 5^^(2) of Eq. ( |1.5| ). /+(0) is calculable in non- 
perturbative QCD. In the flavor SU(3) limit, CVC ensures /+(0)=1, but then ttik = w-tt 
and there are no weak decays. We can write 



/+(0) = l + /2 + /4 + ... 



(3.1) 



In chiral perturbation theory (ChPT) /2 and fi stand for the leading and next-to-leading 
chiral corrections. The Ademollo-Gatto theorem ensures that terms oc {nis — ttIu) are 
absent for vector transitions and /2 = —0.023 is unambiguously predicted in ChPT. The 
calculation of the chiral loop contribution, A(/i) in 



(3.2) 



has been recently completed in ref. |^^, but the full determination of /4 necessitates an 
accurate estimation of the local counter-term /4|^°^(/i), which is 0{p^). Many theoretical 
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Figure 5: Present determinations of /+(0) = (0) |||, |9H^, IT) from lattice QCD and other 

approaches. Hints from ChPt are exploited from all these approaches. 



approaches have been attempted over the years |41], essentially confirming the original 



estimate by Leutwyler and Roos which was obtained in a simple quark model |38]. The 
benefit of these new results, obtained using more sophisticated approaches, lies in the 
fact that we are nowadays in the position to control the systematic uncertainties of our 
calculations while with the quark models this is not possible. To stress the importance of 
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Figure 6: Summary of Jk/ fir estimates All values are from Lattice QCD. In 

recent studies, sea quarks are getting light and data are matched to ChPt fits to determine the 
Low-energy-Constants (LEG). 



the accurate determination of f\, we should remind the reader that the experimental error 
on |I4s|/+(0) is only 0.2%, whereas the spread of theoretical estimates of /+(0) is still at 
the 1% 2% which is unsatisfactory. 

Recent progress in lattice QCD gives us more optimism as far as the prospects of reduc- 



ing the error on /+(0) to well below 1% are concerned [46|. Most of the currently available 
results obtained by using lattice QCD worked with the "heavy pions". Nevertheless, we 
can already see that the lattice QCD results are systematically lower than those obtained 
by the ChPT-inspired models. An important step to resolving this issue has been recently 
made by the UKQCD-RBC collaboration [^9|. Their preliminary result /+(0) = 0.964(5) 
is obtained from the unquenched study with Np = 2 + 1 flavors of the quarks which have 
a good chiral properties on the lattice with finite lattice spacing (so called, domain wall 
quarks), and their pions (> 300 MeV) are much lighter than what is reported in the pre- 
vious lattice QCD studies. Their overall error is estimated to be 0.5%, which is very 
encouraging. It is important to emphasize that they observe the chiral logarithmic cor- 
rections, those which appear in the form factor as /2, which was never observed before 
(most probably because other studies were restrained to heavier pions). One should keep 
in mind, however, that present study is performed at a single value of the lattice spacing 
(i.e. a = 0.12 fm) and in a relatively small lattice box. 

3.2.2 /aV/tt determination 

As we can see in eq. ( |1.6| ), the QCD uncertainty enters with 

Ik/U = 1 + r2 + . . . (3.3) 

In contrast to the semileptonic case discussed above, the Ademollo-Gatto theorem does 
not apply in this case and r2 is not predicted unambiguously in ChPT. Rather one should 
fix the low energy constants from the lattice QCD studies of fx/ fir @1-|45]. Such obtained 
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values are summarized in fig. ^ from whicli we deduce tliat tlie present overall accuracy 
is about 1%. The novelty are the new lattice results with Np = 2 + 1 dynamical quarks 
and pions as light as 280 MeV |^], obtained by using the so-called staggered quarks, 
in which they covered a broad range of lattice spacings (i.e a=0.06 and 0.15 fm) and kept 
sufficiently large physical volumes (i.e. m,r-^ ^ 5.0). It should be stressed, however, that 
the sensitivity of /x/Ztt to the lighter pions is larger than in the computation of /+(0) and 
that the chiral extrapolations are far more demanding in this case. Further improvement 
is expected soon. PACS-CS Collaboration |47] has recently presented preliminary results 
for Np = 2 + 1 clover quarks with pion masses > 200 MeV {m.,^L > 2.9 and a = 0.09 fm). 
For the moment, these simulations are on a single value of lattice volume and finite volume 
effects can be large. 



3.2.3 A test of lattice calculation: the Callan-Treiman relation 



As described in L2 the Callan-Treiman relation fixes the value of scalar FF at t = Akj^ 
(the so-called Callan-Trieman point) to the ratio of the pseudoscalar decay constants fx/ fir 
and the recent parametrization for the scalar FF|p allows this constraint to be exploited. 
The net result is that the ratio (/K//7r)//+(0) can be determined from the measurement 
of the scalar form factor. 



Very recentely KLOE |34] and NA48 [36| have presented results on the scalar form 



factor slopes using the dispersive parematrization. The older results of KTeV and ISTRA+ 
measurement of the scalar form factor slope performed using the 1st order Taylor expansion 
parametrization can be translated in the dispersive parametrization. The results are given 
in Table ^ for all 4 experiments in the case of pole parametrization for the vector form 
factor. The original KLOE and NA48 results are also shown for comparison. 



Experiment 


log(C) direct 


log(C)t 


KTeV 




0.203(15) 


KLOE 


0.207(24) 


0.207(23) 


NA48 


0.144(14) 


0.144(13) 


ISTRA+ 




0.226(13) 



^ Estimated from Ao published. 



Table 9: Experimental results of for log(C) and A+ 



Figure ^ shows the values for /+(0) determined from the scalar form factor slope 
measurements obtained using the Callan-Treiman relation and fx/ fir = 1-189(7). The 
value of /+(0) = 0.964(5) from UKQCD/RBC is also shown. 

As already noticed in ^ the NA48 result is difficult to accommodate, and once compared 
with theory it violates the Fubini-Furlan theorem /+(0) < 1. For this reason the NA48 
result will be excluded when using the Callan-Treiman constraint. 

The average of the experimental results on the EE's with the pole parametrization for 
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the vector FF and the dispersive parematrization for the scalar form factor gives: 



Ai = 0.0256 ± 0.0002 

(3.4) 

A[5 = 0.0149 ± 0.0007 

with correlation coefficient —0.32. The above result is then combined to the lattice deter- 
mination of fK/f-K = 1.189(7) and /+(0) = 0.964(5) using the constraint given from the 
Callan-TVeiman relation. The pole parematerization has been used to describe the vector 
form factor. The results of the combination are given in table ^ 



K 




/+(0) 


fx/ flT 


0.02563(20) 


0.01455(51) 


0.963(44) 


1.1913(61) 


correlation matrix 


1. 


-0.24 


0.11 


-0.13 




1. 


-0.45 
1. 


0.54 
0.27 
1. 



Table 10: Results from the form factor fit. 



where log C = AgAi^^/m^ + 0.0398 ± 0.0041. 

The fit probability is 39% confirming the agree- 
ment between experimental measurements and lat- I 
tice determination. The accuracy of fx / f-w is slightly 
improved. The improvement is better seen in the | 
ratio /+ (0) / {Jk /fn), directly related to the Callan- 
Treiman constraint. This latter improvement is 
very effective in the constraining scalar currents. 

3.3 Test of Cabibbo Universality or CKM ' 
unitarity 

To test CKM unitarity we use the value f+{0)Vus = 
0.21661(46) given in Table |, Ks/Kd/x//. = 
0.27599(59) (see [3T2D , /+(0) = 0.964(5), and 
fK/fn = 1.189(7). The results are: 



Vus = 0.22461 ± 0.00124 (3.5) 
Vus/Vud = 0.23211 ± 0.00145 (3.6) 



ISTRA+ 



UKQCD/RBC 



lavi 

'net 



Kaon WG 



0.9 0.95 1 1.05 

Figure 7: Values for /+(0) determined 
from the scalar form factor slope us- 
ing the Callan-Treiman relation and 
fji/f^ = 1.189(7). The UKQCD/RBC 
result /+(0) = 0.964(5) is also shown. 



These determinations can be used in a fit together with the the recent evaluation of \Vud\ 
from 0+ 0+ nuclear beta decays. The fit gives \Vud\ = 0.97417(26) and = 0.2253(9), 
with x^/iidf = 0.65/1 (42%). The unitarity constraint can also be included, in which case 
the fit gives Vus = sin % = A = 0.2255(7) and xV^df = 0.80/2 (67%). Both resuhs are 
illustrated in Fig. 
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As described in the introduction the test of 
CKM unitarity can be also interpreted as a test 
of universahty of the lepton and quark gauge cou- 
phng: 



CKM 



Gl or \Vud\^ + \Vus\^ + \Vub\^ 



(3.7) 



= 1 and Gp = G^. 

Using the results of the fit we obtain: 

GcKM = (1.16624±0.00039) x 10"^ GeV'^ (3.8) 

In perfect agreement with the value obtained from 
the measurement of the muon lifetime: 

G^ = (1.166371 ±0.000007) x 10"^ GeV-^ (3.9) 



0.230 
>= 



lavi 
■nef^ 

f+(0) = 0.9644(49) 
fK/f.,= 1.189(7) 




with unitarity 



0.225! 



0.970 



0.975 



Figure 8: Results of fits to \Vud\^ \Vus 
and \Vus\/\Vud\- 



The current accuracy of the lepton-quark universality, sets important constraint to the 
model builders of the beyond SM physics scenarios. For example, the presence of a Z' (see 
Fig. ^, left) would affect the universality: 



Gn. = GcKM 



l- 0.007Qe.(Q..-Q..)^^^?^#^ 



(3.10) 





/ s, d \ ^ s, d 

Figure 9: Z' and Higgs exchange. 

In case of Z' in SO(IO) grand unification theories, (QeL = Q/^L = —'^QdL = 1) we 
obtain mz' > 700 GeV at 95% CL, to be compared the one set through the direct collider 



searches, mz' > 720 GeV |2q] . 



3.3.1 Bounds from helicity suppressed processes 

A particularly interesting observable is the ratio of the Vus values obtained from helicity sup- 
pressed processes to that obtained from helicity allowed modes: i?/23 = Vus{K£2) /VusiKi^) . 
According to|45, ^] the value of R123 would be affected by the presence of scalar density or 
extra right-handed currents: 

R123 = 1 + ^^^^^23 (3.11) 

To improve the accuracy of the determination of R123 we use the the values of /+(0) and 



fK/fiT obtained in 3.2.3 ^. In addition, in this scenario both 0^ 0^ nuclear beta decays 



^Ri23 depends only on the ratio ,/'+(0)/(/K//,r 
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and ir^3 are not affected and the unitarity constraint for these modes can be used. The fit 
described in the previous section has been performed assuming unitarity and ahowing for 
two different values of Vus from helicity suppressed(iC — > fiu) and ahowed modes {K — > inv, 
0+ 0^ nuclear beta decays). We obtain: 



Ri23 = 1.00 28 ± 0.0059 



(3.12) 



If the scalar current is due to charged Higgs exchange as shown in Fig. ^, right: 

6Ri23 = (1 - tan^P , n L . o ) (3-13) 



' I + Omtan l3' 

and the measurement of R123 can be used to set bounds on the charged Higgs mass and 
tan p. Figure 10 shows the exlcuded region at 95% CL in the charged Higgs mass-tan /3 
plane. 

The measurement of BR(i? — > tj/)|E9| can be 



also used to set bound on the charged Higgs mass- 
tan P palne. While the B ^ tv can exclude quite 
an extensive region of this plane, there is an uncov- 
ered region in the exclusion corresponding to the 
change of sign of the correction. This region is fully 
covered by the K ^ fiu and the positive correction 
solution for the B ^ tv is fully excluded. 

3.4 Test of Lepton Flavour violation 

3.4.1 Lepton universality and K13 decays 

Search for lepton flavour violation (LFV) in the 
semileptonic decays Ke3 and K^3 is a test of the 
vector current of the weak interaction. It can there- 
fore be compared to LFV tests in r decays, but is 
different to LFV searches in tti2 and K12 decays. 




■ 95% CL from K^fiv lut^^v 
95% CL from B ^tv 
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Figure 10: Exlucluded region in the 
charged Higgs mass-tan /3 plane. The re- 
gion excluded hy B ^ tv is also indi- 
cated. 



The results on the parameter r^e = R 



Exp 



1.0040 ± 0.0044, 



IS 



(3.14) 



in excellent agreement with lepton universality. Furthermore, with a precision of 0.5% the 
test in Ki3 decays has now reached the sensitivity of r decays. The accuracy of this test 
slightely improves using FF's slope values of section 3^2^ 



0.9998 ± 0.0040 



(3.15) 



3.4.2 Lepton universality tests in Kf^, decays 

The ratio Rk = r(i^e2)/r(i^^2) can be precisely calculated within the Standard Model. 
Neglecting radiative corrections, it is given by 



R 



(0) 
K 



m, 
m 



(m 



K 



2 

e 

I (ml: 



m 



2\2 



2.569 X 10-^ 



(3.16) 
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and reflects the strong helicity suppression of the electronic channel. Radiative corrections 
have been computed within the model of vector meson dominance , yielding a corrected 
ratio of 

Rj^ = rP(1 + = 2.569 X 10"^ x (0.9622 ± 0.0004) = (2.477 ± 0.001) x 10"^ 

(3.17) 

Because of the helicity suppression of Ki,2 in the SM, the decay amplitude is a promi- 
nent candidate for possible sizeable contributions from new physics beyond the SM. More- 
over, when normalizing to K^2 decays, it is one of the few kaon decays, for which the 
SM-rate is predicted with very high accuracy. Any significant experimental deviation from 
the prediction would immediately be evidence for new physics. However, this new physics 
would need to violate lepton universality to be visible in the ratio Ke2/K^2- 

Recently it has been pointed out, that in a SUSY framework sizeable violations of 
lepton universality can be expected in K12 decays [pO|]. At tree level, lepton fiavour violating 
terms are forbidden in the MSSM. Loop diagrams, however, should induce lepton flavour 
violating Yukawa couplings as Iv^ to the charged Higgs boson . Making use of 

this Yukawa coupling, the dominant non-SM contribution to Rk modify the ratio to 



pLFV 



TjSM 



1 + 



m 



K 



(3.18) 



The lepton flavour violating term A13 should be of the order of 10~^ — 10~^, as expected 
from neutrino mixing. For moderately large tan /? and Mfj± , SUSY contributions may 
therefore enhance Rk might by up to a few percent. Since the additional term in Eqn. |3.18 
goes with the forth power of the meson mass, no similar effect is expected in Tri2 decays. 

The world average result for gives strong 
constraints for tan f3 and Mfj± (Fig. |ll| (left)). For a 
moderate value of A13 ~ 5 x 10~^, tan /? > 50 is ex- 
cluded for charged Higgs masses up to 1000 GeV/c^ 
at 95% CL. 



4. Conclusions 

Many new precise mesaurements about Ki^ and Ki2 
decays properties have been performed recently al- 
lowing precise tests of the Standard Model to be 
performed. We determine: 




200 
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Charged Higgs mass (GeV) 



1000 



Figure 11: Exclusion limits at 95% CL 
on tan/3 and the charged Higgs mass 
Mh± from \Vus\Ki2/\Vus\K£3 for differ- 
ent values of A 13. 

using /k/ fn = 1.189(7) and /+(0) = 0.964(5) from recent lattice evaluation we obtain: 



/+(0) X Vus = 0.21661(46) 
fx/f. X Vus/Vud = 0.27599(59) 



ud 



0.2246(12) 
0.2321(15) 



(4.1) 
(4.2) 
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in very good agreement with \ Vud\ = 0.97417(26) from 0^ 0^ nuclear beta decays. These 
determinations can be used to evaluate the Fermi constant from the quark sector: 

GcKM = (1.16624 ± 0.00039) x 10"^ GeV-^ 

to be compared with that obtained from the muone hfetime: 

= (1.166371 ±0.000007) x 10"^ GeV"^ 

The improved accuracy of these measurements sets non-trivial constraints for physics be- 
yond the Standard Model. In particular, the comparison of the \ Vus\ determinations from 
helicity suppressed and allowed processes sets a competitive lower bound on the charged 
higgs mass as a function of tan f3. 

Very recentely the test of lepton flavour violation in helicity suppressed processes 
through the observable Rk = T{Kf.2)/^{Kn2) as been improved significantly. The new 
world average 

Rk = BR{Ke2)/BR{K^2) = (2.45 ± 0.11) x 10"^ 
gives the best constraint for the e — r LFV coupling of MSSM. 
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